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We study interacting one dimensional (ID) quantum lattice gases with inte- 
grable impurities. These model Hamiltonians can be derived using the quan- 
tum inverse scattering method for inhomogeneous models and are by con- 
struction integrable. Absence of backscattering at the impurities is shown to 
be the characteristic feature of these disordered systems. The value of the 
effective carrier charge and the Sutherland-Shastry relation are derived for 
the half-filled XXX model and are shown to be independent of the impurity 
concentration and strength. For the half-filled XXZ model we show that there 
is no enhancement of the persistent currents for repulsive interactions. For 
attractive interactions we identify a crossover regime beyond which enhance- 
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Abstract 



merit of the currents is observed. 
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Much work has been done recently towards a deeper understanding of the effect of impu- 
rities in Fermi- and Luttinger hquids. Conformal field theories with boundaries have been 
used to study the multichannel Kondo problem [|^. Bosonization techniques and renormal- 
ization group methods have given new insights into the problem of potential scatterers in 
Luttinger liquids [0-^. Also, the observed discrepancy between the experimental value of the 
persistent current in moderately disordered ensembles of quasi one-dimensional (ID) rings 
P and the theoretical predictions based on calculations in a disordered but non-interacting 
electron gas P|, have generated much interest in the interplay between interactions and 
impurities [[^-[TT||. 



Besides the spin impurity model of Andrei and Johannesson [0, no other exactly solvable 
interacting quantum many-body problems with impurities have been studied. This is in 
retrospect somewhat surprising as the existence of models with site-dependent defects has 
been mentioned throughout the literature on the method of quantum inverse scattering 



(QISM) [p!3| -[r^. Recently, Bares ||I8|] and, independently, Schmitteckert, Schwab and Eckern 
(SSE) |TP| have used these remarks and explicitly investigated a tJ model and an XXZ model 
with integrable defects, respectively. 

In the XXZ model, SSE find the following surprising properties: (i) The energy spectrum 
is independent of the spatial distribution of impurities; (ii) there is no localisation of the 
ground state wave function. In contrast measurements of the transport properties in meso- 



scopic systems show a very sensitive dependence on the impurity configuration ||2T]]. Further, 
it is well understood that in the thermodynamic limit in ID a generic impurity drives the 
system from a metallic to an insulating phase even in the presence of weak interactions 

iiil- 

In the present work, we shall first introduce the (antiferromagnetic) XXX Heisenberg 
model as a simple example of a system with integrable impurities. We will then show that a 
complete physical understanding of the anomalous properties of this class of models can be 
got by studying the properties of a single particle scattering off a single impurity. The system 
is then threaded with an external flux $. For the half-filled band, we establish the charge 



of the effective carriers to be —q/^. The Sutherland-Shastry relation p3|j2^ Dx = 1/2% 
still holds in the presence of impurities. For the XXZ model with integrable impurities at 
half filling the charge stiffness D shows enhancement of the persistent currents for attractive 
interactions in agreement with studies of XXZ models with potential impurities M. 



We construct our model from the R and L„ matrices of the Heisenberg XXX model 
i.e. R{X) = {XP + i)/{X + i) and L„(A) = A + | ELi o"" ® <• P is a 4 x 4 permutation 
matrix and {cr'*} are the Pauli matrices. These operators satisfy the local Yang-Baxter (YB) 
equation, i.e. R{X—X')[Ln{X)®Ln{X')] = [Ln{X')^Ln{X)]R{X—X'). The main idea that allows 
for the construction of models with integrable impurities is that the YB equation continues 
to be satisfied under an arbitrary local shift in the spectral parameter Ln{X) — > L„(A + Un). 
The transfer matrix T with local shifts {z/i, ■ ■ ■ , uj^} is given as T(A, {i^n}) = Iln=i LniX + Un). 
The Hamiltonian H is then given as usual by the logarithmic derivative of the transfer matrix 
T(A, {i^n}) evaluated at the special point Aq = Note that this construction is different 



from the one used by Andrei and Johannesson |T2| and does indeed yield inequivalent, albeit 
similar, Bethe Ansatz (BA) equations. 

As is customary, we now Jordan- Wigner transform the spin model into a lattice gas of 
spinless fermions and impose periodic boundary conditions. The resulting Hamiltonian on 
a ring of size with M spinless fermions and a single impurity at site m with strength u is 
given as i^" = Hxxx + Hm,imp, where Hxxx = - Y.f=i{c]+iCj + h.c.) + 2 Y.f=i PjPj+i + 2M is 
the usual Hamiltonian of the clean XXX model in units of the hopping energy. The operators 
c] and Cj create and annihilate fermions on site j and pj = c^jCj is the density operator. The 
impurity interaction at site m is Hm,imp = Hm,t + H^^p + Hmj with 

Hm,t = — (m + iw)[cl^Cm.-l + c'm+l'^ni + <^rn~l(^m+l\ 

+ h.c, (la) 

H^^p = -2u[p 

m Pm (Pm-1 + Pm+l) + Pm-lPm+l], (lb) 
Hm,j ~ 2lt(c^_(_]^Cm— 1 ~l~ C^m—l^rn+l) Pm 

-\- 2iw(cly^Cm—l (^rn-l^rn) Pnri+l 



where the couphng constants as a function of u are u = —v^j (1 + z^^) and w = —v /{I + v"^)- 
In Eq. (1), we have three types of terms: (i) Hra,t'- Nearest and next nearest neighbour 
hopping matrix elements are modified involving the sites m and mil. (ii) Hm,p'- The onsite 
potential at site m and both nearest and next nearest neighbour density-density interaction 
terms are induced, (iii) Hm,j'- This term corresponds to a current density interaction. Thus 
we see that a local translation of the spectral parameter A A + z/^ effectively involves the 
site m and its two neighbors mil and in addition involves next nearest neighbour terms 
involving the sites m — 1 and mil. We note that the impurity terms involve qualitatively 



the same processes as in the tJ model with integrable impurity ||18 . 

The Hamiltonian in the presence of more than one impurity can now easily be con- 
structed: In case the impurities are well separated (|m — m'| > 1) the Hamiltonian reduces 
to a sum over the isolated impurities i.e. Y^a=i Hm{a),imp for F such well-separated impurities 
at positions m(a), a = 1, . . . , F. In case two impurities are on neighboring sites, say m and 
mil, the Hamiltonian will acquire similar terms involving the sites m — 1 to m i 2. As we 
will show later, this does not change any of the physics of these systems. 

To understand the properties of integrable impurities, we first study the scattering of a 
single particle by a single impurity at site m. The terms in the Hamiltonian involving only a 
single particle then are = - J2i{clci+i+h.c.)+Ht-2upm+'2. Let \E'i„(n) = e*''"ii?(/c)e"*''"' 
for n < m be the incoming and "^outiiT-) = T{k)e^^^ for n > m he the outgoing wave 
functions and R{k) and T{k) the reflection and transmission amplitudes, respectively. Then 
with e{k) = — 2(cosA; i 1) the usual nearest-neighbor one-particle dispersion, we find that 
(i) R{k) = for all momenta k. Thus there is no backs cattering of a particle due to the 
impurity. This is consistent with the understanding that the number of scattering channels 
in a BA-solvable Hamiltonian is conserved. In contrast, a generic impurity would give rise 
to a reflected and a transmitted wave thus doubling the number of scattering channels. 
We remark that the parametrisation of u and w given in Eq. (1) is unique, i.e. it is the 
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only such parametrisation that will lead to R{k) = for all k and u. (ii) Unitarity of the 
scattering matrix requires |Tp + = 1 and we consequently have \T{k)\'^ = 1. Writing 
T{k) = e*'^(^''^\ we then find the phase-shift of the particle wave function due to the presence 
of the impurity to be 6{k, v) = —k + 2 arctan(tan(A;/2) — 2z/). 

To study the problem with more than one particle, we use the BA equations in the 



presence of the impurities which are given for the XXX model by QISM as [|I4 

M F 

Nkj = 2nlj - 9{kj, kk) - ^(kj, Umia)), (2) 

k=l a=l 

where 6{k) = 2 arctan[(tan(fc/2) — ta.n{k' /2))/2] is the two-body phase shift. Note that 
although the total momentum P is no longer a good quantum number in the presence of 
impurities, the individual BA pseudo-momenta k still represent the eigenvalues of the BA 
equations and are thus good quantum numbers. 

We now make the important observation that the only change in the BA equations w.r.t. 
the clean case is the additional phase shift 6{k, v) which is identical to the forward scattering 
phase shift acquired by a single particle scattering off the impurity. This is true although 
the BA equations include many-body interactions and the interactions Hm(a),p and Hm{a),j 
induced by the impurities. Thus it is also irrelevant, if the impurities are well separated 
or not. For all configurations of the impurity the net effect will be the same phase shift 
6{k, for each impurity of strength z/^. Further, Eq. (Q) is insensitive to the spatial 
distribution of the F impurities. This is a consequence of the absence of backscattering. 
Quantum interference of multiply backscattered waves would in general lead to a non trivial 
dependence of the energy spectrum on the impurity configuration. Note that, although the 
energy spectrum is insensitive to the impurity configuration the wave function is no longer 
given by the simple coordinate BA form '^{xi, . . . ,xm) = ^(P) exp(i Y^fLi Xjkpj) with 
P a permutation of the M momenta kj. Presumably, one might calculate the wave function 
from the reference state of the QISM, which then by construction includes a site dependence. 

We believe that these results are generic for all such integrable impurities constructed 
by QISM. This is implied by the very existence of the BA equations. Besides the XXX 



model, the tJ model of Bares [T^ and the XXZ ring of SSE pO[, this should therefore also 
apply to, e.g., the Hubbard model with integrable impurities constructed from two coupled 
inhomogeneous six- vertex models |T^ . 



To study the transport properties of these integrable impurity models, we thread the 
system with an external flux $, which accelerates all the particles around the ring. The 
variation of the energy spectrum as a function of flux is made nontrivial for an interacting 
quantum system only in the absence of Galilean invariance. In the present problem this 
invariance is destroyed both due to the presence of the lattice and the existence of the 
impurities. 

The ground state energy -Eo($) for a given flux $ is periodic in $ with period 2tt. 
However, the adiabatic variation of the energy £'($) of the state that begins as the ground 
state at $ = may have a periodicity greater than 27r. The clean {v = 0) XXX model 
has been studied previously by Sutherland and Shastry ||2^. They find a periodicity of 4% 



implying carriers with charge —q/2 in terms of the charge q of the fundamental carriers. 
The second quantity of interest is the stiffness D = {N/2)d'^EQ{^)/d^'^\,^=o. Disorder will 
in general destroy the phase sensitivity of the ground state and therefore lead to a non trivial 
renormalization of the stiffness. 

Let us consider the half-filled case = 2M. In order to study the periodicity of the 
ground state, we follow the particle with the maximum momentum fcjv/ characterised by the 
ground state quantum number Im = {M — l)/2 as a function of $. At the special point 
f^M = TT, the BA Eq. for kM with 6{ki, tt) = tt for i < M and 6{tt, u) = for all finite u 
gives the flux to be 2tt. Conversely, we notice that as long as I'I'I < 27t, all k's stay within 
\k\ < TT until at $ = 2tt, the maximum momentum takes the value k^j = vr. Hence, it follows 
that for 27r < $ < 47r the k's are the same as in the case |$| < 27r with the states relabelled 
as kM~i > > ■ ■ ■ > ki > kjif modulo 2tt. We next substitute the values kM = and 

$ = 27r into the remaining M — 1 equations. These reduce to 

F 

Nk^ = 27r/; - J2 ^ik,, h) - J2 S{k„ z/„.(„)). (3) 

l^M a=l 



We immediately notice that these equations with the quantum numbers Jj = Ij + ^ are just 
the BA equations for the ground state of M — 1 particles in the absence of flux. The energy of 
the ground state of the M particle system E{^) at flux $ = 27r is thus related to the ground 
state of M - 1 particles at $ = 0, i.e. ^M(27r) = ^o,m-i(0) - 2(1 + cos(A;m)) = ^o,a/-i(0). 
Therefore 

AEm{^) = Em{27t) - Em{0) = Att^D/N 

= Eo,M-iiO) - Eo,m(0) = 2nx-'/N, (4) 

with X the susceptibility. This remarkable property of the BA Eq. (||) relates the stiffness 
and the susceptibility through Dx = l/27r even in the presence of impurities. This is the 



Sutherland-Shastry relation derived previously in the clean XXZ Hubbard ||2^ and SC 
model 1^ . 



We next turn to explicitly calculating the stiffness D of a system with integrable impu- 
rities. For convenience, we will use the XXZ spinless fermion model with impurities where 
we can easily vary the interaction strength A = cos{2ri). The clean system is given by the 
Hamiltonian, Hxxz = — Yl^=i{c]+iCj + h.c.) — 2A SjLi PjPj+i- The stiffness is calculated by 
studying the leading order flnite size corrections to the thermodynamic value of the ground 
state energy in the presence of a flux with the help of the Wiener-Hopf technique ||2^. For 



the XXZ model, this was already done by SSE for impurities chosen symmetric w.r.t. the 
origin and of equal strength u. They get for the stiffness 

r,/ X ^ 7r sm{2r])/{7]{7i -27])) 
"^^"l^^) 8 1-/ + /cosh[7rz//(7r-2r7)]" 

where / = E/N is the impurity density. By the same arguments as above, Dx = 1/2tt also 
in the XXZ model and we have the susceptibility in the presence of impurities. Using hydro- 
dynamical arguments, we further have Dx~^ = |2^, where v is the Fermi velocity. Thus 



we see that the impurities simply renormalize the values of the stiffness, inverse susceptibil- 
ity and the Fermi velocity of the clean system by the factor (1 — / + / cosh[7rz//(7r — 27/)])^^. 
This then further implies that the critical exponents remain unchanged and that the corre- 
lation functions are modifled by the renormalization of v only. We remark here that these 
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calculations can also be done by choosing a more general, say Gaussian, distribution for the 
impurity strengths. However, the resulting expressions are lengthy and will not be included 
in this letter. 

A finite stiffness also implies that the ground state is extended. The stiffness in the 
z/ 7^ case as seen from Eq. (^) has a quadratic dependence on the strength of the impurity 
and only in the limit of very large disorder does one get an exponential dependence. This is 
contrary to the understanding that due to coherent back scattering in ID a generic impurity 
would lead to localisation for arbitrarily weak disorder. However this result is consistent 
with the observation that these are "special" impurities with zero reflection. 

The stiffness provides an operational deflnition for the persistent currents as J = D$. 
In Fig. |l|, we plot D of the half-fllled XXZ model with impurity density / = 0.001 in 
the attractive and the repulsive interaction regimes as a function of impurity strength u. 
In agreement with previous studies of the XXZ model in the presence of scalar potential 
scatterers ||^, we flnd that there is no enhancement in the persistent currents for repulsive 
interaction (A < i.e. 7r/4 < t] < 7r/2). For attractive interactions (A > i.e. < ?7 < 
71/4), we observe such a crossover into a regime where the currents for non-zero attraction r) 
are larger than the currents of the non-interacting case at the same strength of the impurities 
0. Thus, although there is no backscattering in the present case of integrable impurities, 
the situation regarding enhancement of persistent currents is qualitatively the same as for 
potential scatterers. 

Away from half-flUing, the situation is more complicated. Here we only remark that the 



degeneracy of levels observed in the clean XXZ for 1/r flllings with r rational is lifted. 
Details will be published elsewhere. 

In conclusion, we have shown that an interacting quantum many-body system with 
impurities, that can be constructed by the QISM, is characterised by an absence of reflection 
at the impurities. As a direct consequence of the absence of back scattering, the ground state 
of the interacting system remains extended even in the presence of these special impurities. 
These systems are therefore unique in that there is no localisation in the presence of disorder 



in ID. The value of the effective carrier charge — g/2 and the Sutherland- Shastry relation 
Dx = l/27r are derived for the half- filled XXX and XXZ models and are shown to be 
independent of the impurity strength. For the half-filled XXZ model, we show enhancement 
of the persistent currents only for attractive interactions. 
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respectively. 



10 



REFERENCES 



* Permanent address: Institut fiir Physik, TU Chemnitz- Zwickau, D-09107 Chemnitz, 
Germany. 

[1] I. Affleck, in Springer Series in Solid-State Sciences 118, ed. by A. Okiji and N. 

Kawakami, (Springer Verlag 1994). 

[2] T. Giamarchi and H. Schulz, Europhys. Lett. 3, 1287 (1987). 

[3] T. Giamarchi and H. Schulz, Phys. Rev. B 37, 325 (1988). 

[4] C. L. Kane and M. P. A. Fisher, Phys. Rev. Lett. 46, 15233 (1992). 

[5] L. P. Levy G. Dolan, J. Dunsmuir and H. Bouchiat, Phys. Rev. Lett. 64, 2074 (1990). 
V. Chandrasekhar, R. A. Webb, M. J. Brady, M. B. Ketchen, W. J. Galager and A. 
Kleinsasser, Phys. Rev. Lett. 67, 3578 (1991). D. Mailly, C. Chapeher and A. Benoit, 
Phys. Rev. Lett. 70, 2020 (1993). 

[6] B. Altshuler, Y. Gefen and Y. Imry, Phys. Rev. Lett. 66, 88 (1991). 

[7] R. A. Romer and A. Punnoose, to be pubhshed in Phys. Rev. B (1995). 

[8] T. Giamarchi and B. S. Shastry, Phys. Rev. B 51, 10915 (1995). 

[9] G. Bouzerar, D. Poilblanc and G. Montambaux, Phys. Rev. B 49, 8258 (1993); M. 
Ramin, B. Reulet and H. Bouchiat, Phys. Rev. B 51, 5582 (1994). 

[10] M. Kamal, Z. H. Musslimani and A. Auerbach, preprint (1995). 

[11] H. Mori and M. Hamada, preprint (1995). 

[12] N. Andrei and H. Johannesson, Phys. Lett. A 100, 108 (1984); K.-J.-B. Lee and P. 
Schlottmann, Phys. Rev. B 37, 379 (1988); P. Schlottmann, J. Phys. CM 3, 6617 (1991); 
E. S. S0rensen, S. Eggert and L Afffeck, J. Phys. A 26, 6757 (1993). 

[13] R. J. Baxter, Exactly solved models in statistical mechanics, (Academic Press, London, 

11 



1982). 

[14] L. A. Takhtajan, in Vol 242 Exactly Solvable Problems in Condensed Matter and Rel- 
ativistic Field Theory, edited by B. S. Shastry, V. Singh and S. S. Jha (Springer, New- 
York, 1985). 

[15] B. S. Shastry, J. of Stat. Phys. 50, 57 (1988). 
[16] E. K. Sklyanin, preprint HU-TFT-91-51 (1991). 

[17] V. E. Korepin, N. M. Bogohubov and A. G. Izergin, Quantum Inverse Scattering Method 
and Correlation Functions, (Cambridge, New York, 1993). 

[18] P. A. Bares, preprint (1994). 

[19] P. Schmitteckert, P. Schwab and U. Eckern, Europhys. Lett. 30, 543 (1995). 

[20] Note that for example in the XXZ model of SSE the impurity phase shift 5 is not only 
a function of k and u but also of the interaction strength A. 

[21] B. L. Altshuler, P. A. Lec and R. A. Webb, (eds.), Mesoscopic Phenomena in Solids, 
(North- Holland, New York, 1991). 

[22] N. F. Mott and W. D. Twose, Adv. Phys. 10, 107 (1961) R. E. Borland, Proc. R. Soc. 
London Ser. A 274, 529 (1963). 

[23] B. Sutherland and B. S. Shastry, Phys. Rev. Lett. 65, 1833 (1990). 

[24] B. S. Shastry and B. Sutherland, Phys. Rev. Lett. 65, 243 (1990). 

[25] R. A. Romer and B. Sutherland, Phys. Rev. B 50, 15389 (1994). 

[26] F. Woynarovich and H. P. Eckle, J. Phys. A. 20, L97 (1987); C. J. Hamer, G. R. W. 
Quispel, and M. T. Batchelor, J. Phys. A 20, 5677 (1987). 

[27] B. Sutherland, unpublished; R. A. Romer and B. Sutherland, Phys. Lett. A 190, 295 
(1994); B. Sutherland, R. A. Romer and B. S. Shastry Phys. Rev. Lett. 73, 2154 (1994). 



12 



[28] R. A. Romer, H.-R Eckle and B. Sutherland, Rhys. Rev. B 52, 1656 (1995). 



13 



FIGURES 

FIG. 1. Stiffness of the half-filled XXZ model with impurity density / = 0.001 for various 

interaction strengths A as a function of impurity strength u. Note that increasing repulsion and 
impurity strength will always reduce D whereas there is a crossover for all values of the attraction 
as we increase the impurity strength. 
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